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Abstract. The acceleration of electrons at 1-10 keV energies is the cause of the polar aurora 
displays, and an important factor of magnetic energy transfer from the solar wind to the Earth. 
Two main families of acceleration processes are observed: those based on coherent quasi-static 
structures called double layers, and those based of the propagation of Alfven Waves (AW). This 
paper is a review of the Alfvenic acceleration processes, and of their role in the global dynamics of 
the auroral zone. 
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INTRODUCTION 

The ionosphere is like a screen on which the electrons coming from the magnetosphere 
are projected, giving rise, when they carry enough energy, to visual displays, the polar 
auroras. Among them, dijfuse auroras are very common and hardly noticeable. Discrete 
auroras can be seen in visible light, in UV, and even in X rays. Some of them are on the 
day side, therefore they are rarely seen (except in winter). Some others are very bright 
and appear mainly on the night-side; their appearance is correlated with rather sudden 
changes in the magnetic configuration of the magnetosphere, for instance substorms and 
pseudo-breakups [1]. 

As the other auroras, discrete auroras are visible at altitudes comprised between 80 
and 400 km, in two ovals around the northern and the southern magnetic poles, for 
magnetic latitudes in the range 60-80°. The auroras are visible in the ionosphere, and 
the radiative processes that subtend them are well understood. We know, since the first 
rocket flights in the ionosphere that they are caused by fluxes of electrons with an energy 
of the order of a few keV. This energy per particle is much larger than the thermal energy 
of the ionospheric particles (~ 1 eV) and of the solar wind (~ 100 eV). 

The present study is mainly concerned with the bright discrete auroras that can be 
seen on the night side. Our purpose is to introduce a series of theoretical works about 
the cause of the acceleration of the electrons, that we now call "auroral acceleration" . 
Not all the mechanisms are reviewed. We focus our study on the role that Alfven waves 
(AW) can play in the auroral acceleration. First, a brief overview of the energetics of 
the general magneto spheric processes is given. The importance of the electric current 
systems is emphasized, as an energy source for the auroral acceleration. The next section 




FIGURE 1. Left: Schematic view of the magnetosphere, the main currents and the auroral zone, 
evidenced on this figure through the auroral currents. Right: Solar wind/magnetosphere energy flowchart 
Reprinted with permission from [2]. Copyright [1984], Springer. 



addresses the very important role of quasi-static electrostatic structures in the generation 
of the discrete auroral arcs. Then, the role of the Alfven waves on auroral acceleration 
is considered, paying attention to the frequency range. 



Current systems in the magnetotail 

Contribution of auroras to magnetospheric energy transfers 

The Earth environment receives a continuous flux of energy from the solar wind. 
Stern [2] made a review of the various fluxes of energy in the magnetosphere, from 
the solar wind down to the ionosphere, and evaluated the overall flux of the solar wind 
flowing through the magnetosphere section as 1.3 x 10^^ W from which about 1 percent 
is captured by the magnetosphere. Among the various magnetospheric processes, the 
most energetic is the transverse electric current across the magnetotail, requiring a power 
of 3 X 10^^ W (estimated up to a distance of 40 Earth radii). A large part of the energy 
is accumulated into the magnetotail when the magnetosphere is quiet, while a flux of 
1.3 X 10^^ W is continuously carried by the plasma towards the Earth, across the inner 
boundary of the magnetotail current layer. During substorms, a more substantial fraction 
of it (3 X 10^^ W) is released toward the Earth. This energy corresponds to the liberation 
of the magnetic energy accumulated into the magnetotail during the quiet times, and it 
corresponds also to a full transfer of the energy that transits through the magnetotail 
current. 



Most of the energy directed from the magnetotail towards the Earth is conveyed 
through the field aligned currents (FAC), also called Birkeland currents, which are elec- 
tric currents in the auroral zone propagating along the magnetic field lines. According to 
Stem, the energy flux along the FAC is 1.4 x 10^^ W during quiet times, and 2.7 x 10^^ 
W during substorms. 

He considers 4 processes connecting the FAC and the ionosphere 

1. Some of the FAC are connected to the ring current and can provide energy to it. 
This energy is lost for the ionosphere. (The ring current is a vast electric current 
flowing in the magnetosphere at a distance ~ 8 Earth radii all around the Earth.) 

2. In the auroral zone, some plasma is accelerated towards the Earth. The acceleration 
results form parallel electric fields These processes are very sensitive to the 
magnetospheric activity. Auroral plasma acceleration involves energy fluxes in the 
range 1 — 5 x 10^*^ W, the highest values being reached during substorms. 

3. The same parallel electric field can accelerate a part of the plasma out of the iono- 
sphere. Stern argues that the involved flux of energy are negligible in comparison to 
the others. (Some recent observations and theories put this assertion into question 
as will be shown later in this review.) 

4. The FAC reach the ionosphere. Then, the current system is closed in the ionosphere, 
mainly with horizontal currents propagating in a medium of finite (Pedersen and 
Hall) conductivities. Then, 2 — 16 x 10^° W come into Joule dissipation. 

The sum of all these energy fluxes is, again, comparable to those transiting through 
the internal layer of the magnetotail current sheet. Therefore, the auroral region is the 
region where most of the energy flowing Earthward from the magnetotail is dissipated. 

Auroral acceleration in a finite box 

The study of auroral acceleration is therefore connected to the transit of energy from 
the magnetotail towards the ionosphere. In Fig. 1, this transition is represented by a box 
with an input of field aligned current, and outputs including plasma acceleration and 
auroral kilometric radiation. Within a theoretical approach, it is also necessary to define 
a box, but in the configuration space, with an upward boundary on the magnetospheric 
side, and a downward boundary on the ionosphere, or slightly above, on which boundary 
conditions are to be prescribed. As with electric circuits, we are concerned with a (field 
aligned) electric current and an (accelerating) electric field. With electric circuits, the 
boundary conditions can be set equivalently as a voltage generator and an impedance in 
series (Thevenin circuit), or a current generator mounted in parallel with an impedance 
(Norton circuit). The two systems are equivalent from a computational point of view, but 
the physics behind them is not the same. In the auroral zone, a "normal conductivity" 

is defined by J|| = C7||E|| and C7|| ~ > oo. Because ay oo, it is easy to create a 

current, and difficult to set an electric field. To prove acceleration, one needs to show the 
rise of parallel electric field and a kind of anomalous conductivity. Thus, a forced current 

is to be considered and the question is how it could trigger an accelerating electric field. 
(Setting this current imposes strict conditions on the magnetic field [3].) 



Defining a resistivity may seem paradoxical, because the resistivity that triggers the 
electric field, and therefore the acceleration, tends to be opposed to particle motion. The 
paradox can be solved if we consider that the resistivity imposes a deceleration of most of 
the plasma and the efficient acceleration of a minority of particles that carry the (forced) 
parallel current. This way of resolving the paradox (that is confirmed by observations) 
requires to go beyond the mere distinction between mean and thermal velocities. This is 
why the study of auroral acceleration cannot be based on a fluid theory, but on a kinetic 
(or hybrid) description of the plasma that takes into account minor components of the 
electron distribution function. 

The first example of such a kinetic process was proposed by Knight [4], it is based 
on the increase of the magnetic field amplitude as particles approach the ionosphere. 
Then, the mirror effect comes into action: the particles with a finite pitch angle bounce 
back. Therefore, they cease to carry the field aligned current. If the current is forced, an 
electric field is created : it counteracts the mirror effect for a minority of particles (those 
with low pitch angle) that are accelerated to maintain the current at low altitude. 

Another effect that would favour parallel electric fields is the rarefaction of the 
electron density : it there are less electrons to carry the current, as the current is forced, 
some electrons must be accelerated [5]. 



Time independent j\\ : quasi-static acceleration region 

Knight theory involves two boundaries : the ionosphere and the magnetosphere. On 
the former (/ index) the plasma is cold and the magnetic field high. On the latter (M 
index), the plasma is hot and the magnetic field weaker. The relation between the electric 
potential drop A<i> and j\\ (obtained from the Vlasov-Poisson equation) reads [4] 

for 1 « -— and -— « 7 , ~ -cNm ^ (1) 
Ti Tm Bm " yiTimeJ Tm 

Numerically, the current 7|| needed for the observed 1-10 keV potential drop is larger 
than the measured values by one order of magnitude. Therefore, the convergence of the 
magnetic field lines is favourable to the downward acceleration of charged particles, but 
not sufficient for the observed auroral electron energies. 

Actually, Knight theory does not deal with the reaction of the plasma on the electric 
field. Experimental [6] and theoretical [7, 8, 9, 10] works showed that when an electric 
potential drop appears in a plasma, it tends to concentrate over a region of small extent (a 
few Debye lengths). It forms a coherent electrostatic structure, called double layer (DL) 
which is said to be strong when the potential drop exceeds the plasma thermal energy, 
and weak otherwise. 

In space, the analysis of accelerated particles showed that strong and weak DL exist 
in the auroral zone. Strong DL, for instance, are located at a few thousands of km above 
the ionosphere [11, 12]. The analysis of the electric field reveals their multi-dimensional 
structure, with V-shaped equipotential lines, showing that the associated electric field is 
not systematically parallel to the field aligned current. A sketch of a strong DL in an 
upward FAC is shown on Fig. 2. The continuous lines are electric equipotentials and 



the dotted lines are the magnetic field lines. In the DL, their gradient, and therefore the 
electric field are aligned with the magnetic field. This is the region of acceleration. Other 
details on the figure show various phenomenon that are triggered by the accelerated 
particles, like plasma turbulence, auroral kilometric radiation (AKR) and other waves, 
and plasma cavities. These structures are multidimensional, but DL with non-aligned 
electric fields can be already described with a ID model [13]. A few theories take 
into account the difference of magnetospheric and ionospheric temperatures [14, 15], 
emphasize the role of mass transfer [16], take into account the vertical density gradients 
[17] or the convergence of the magnetic field lines and the effect of various populations 
of particles [18, 19]. Many of these models impose an electric potential drop across the 
box, instead of a forced current. The simulations do not really explain the cause of the 
accelerating electric field, but show how it evolves under the reaction of the accelerated 
plasma. 

Such structures appear both in regions of upward currents with precipitated electrons 
and bright auroras [20, 21, 22] and in the adjacent downward current regions [23, 24, 
20, 25], with precipitated ions and no auroras (termed black auroras by contrast with 
the neighbouring luminous structures). The acceleration structures in downward current 
regions were less expected (by Stem for instance [2]) because Knight effect is not 
favourable to electron acceleration in the upward direction. 

The duration of a strong DL in an downward current during a substorm could be 
estimated with multi-spacecraft observations [20]: the structure grew in size during ~ 
100 s, then dissipated in a few min. This is the time required to evacuate the accelerated 
electrons from the ionosphere. It is also the typical type of evolution of the magnetotail 
current layer during substorms. Similar times of evolution have been observed for strong 
DL in downward current regions. 

Although a double layer is a stationary solution of the Vlasov and Poisson equations, 
stationarity does not imply stability [26, 27]. Newman et al. have shown for instance that 
the interaction of the accelerated plasma with a cold ambient plasma can trigger a Bune- 
man instability that disrupts the DL [28]. The double layer reforms later, but in a more 
turbulent regime [29], and can be destroyed again. Thus, the life duration of a DL does 
not depend directly on those of the field aligned current, but rather on its interaction with 
the ambient plasma. Therefore, we can consider that the quasi-electrostatic acceleration 
structures are associated to a basically direct current i.e. with very low time-dependency. 
This low time-dependency has little influence on the evolution of the quasi- static accel- 
eration structures. 

What is the importance of strong double layers in auroral acceleration ? In the upward 

current regions, large parallel electric fields have been seen at the boundary between 
the ionosphere-dominated plasma and the auroral depleted plasma, and they contained 
roughly 10-50% of the potential [5, 22]. Observations in other contexts suggested that 
they might be dominant. Nevertheless, the DL cannot explain completely auroral accel- 
eration. Newell et al. have made a statistical analysis over 11 years of data [30]. They 
show that the most significant input of energy per precipitated particles is (surprizingly) 
due to the diffuse aurora, that are not especially concerned by the present review. Then, 
comes into play the contribution of the quasi-mono-energetic beams of particles that are 
supposed to come from acceleration by strong double layers (10 % of the energy flux 
during quiet conditions to 15 % during moderately active periods). Then, broadband ac- 
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FIGURE 2. A strong double layer region: sharp potential drops and V structures. Reprinted with 
permission from [19]. Copyright [2002], American Institute of Physics. 



celeration associated to another type of process, that is only 6 % during quiet periods, 
rises to 13 % (almost as much as DL acceleration) during active periods. The broadband 
acceleration by waves represents also to 28 % of the number flux of particles, dominated 
only by the diffuse auroras (48 %). 

AURORAL ACCELERATION BY ALFVEN WAVES 

Alfven waves as 7|| and energy carriers 

If strong double layers associated with quasi-static FAC do not explain all the flux of 
accelerated particles, could time varying FAC induce other acceleration processes ? 
According to the linear theory of MHD Alfven waves, the parallel electric current 

oscillates with an amplitude 7|| = = ^^^J- ^^^d the Poynting flux is 5 1 = —^f- 

Therefore Alfven waves are able to carry a parallel electric field, and an energy flux. 
Several measurements in the auroral zone identify Alfven waves with a Poynting flux 
large enough to allow for the acceleration of the electrons, even for an acceleration 
process with an energetic efficiency of only a few percents [31, 32, 33]. 

For a shear Alfven wave of frequency = > the linear MHD polarization in a 

homogeneous medium is Ei = (— 5oviy, 0, 0), where Bq = (0, 0,5o) ^iid k = {kx-,^-, ^z) = 
(^j^,0,^l|). For magnetosonic waves it is Ei = (0,5o^ix)0)- In all cases, E-^ = E^^ = 0. 



Therefore, a purely MHD wave in a homogeneous medium cannot generate a parallel 
electric field. Therefore, to see how AW can create a parallel electric field, it is necessary 
to consider a non homogeneous medium, or to go beyond MHD, or both, as we will see 
in the following sections. 

Magnetic fluctuations induced by dtj\\. Alfven waves ? 

Time varying FAC were inferred well before the space age, thanks to the observations 
of the magnetic micro-pulsations. These are magnetic fluctuations registered from the 
ground. They have now been classified, relatively to their period and to the shape of the 
wave envelop. Pil (>100s) and Pi2 (<40s) are irregular micropulsations [34]. They are 
correlated with magnetic fluctuations in space [35], with substorm onsets (start of bright 
and active auroras) [36, 37] , and with current disruption in the magnetotail [38]. From 
a theoretical point of view, the irregular Pil are seen as Alfven wave (AW) field line 
resonance (FLR). Their frequencies correspond quite precisely to the values expected 
for the oscillation of a whole terrestrial dipolar field line. The irregular Pi2 are seen 
as FLR with nodes at low altitude [39, 40]. The bursts of pulsations PiB (0.2-1 Hz) 
correspond to AW trapped at low altitude [41, 42, 43]. The other pulsations (PCI, pearls 
etc) are associated to other kinds of waves or instabilities, without direct link to AW. 

Can these Alfven waves accelerate electrons ? If it is the case, Alfven wave accel- 
eration processes are possibly based on various time-scales: (i) long wavelength AW : 
magnetic field line resonance (FLR). Period > 100s; (ii) meso-scale AW, period <40s, 
often 1 s. They are the resonance of a limited part of a magnetic field line. They are 
associated to the Ionospheric Alfvenic Resonator (lAR); (iii) Bursty wave packets. Du- 
ration < 1 second; (iv) Not seen from the ground : Series of wave packets named Solitary 
kinetic Alfven Waves (SKAW), with a duration of about 10 ms by wave packet. 



In-situ observations of accelerated particles 

The in-situ observations of accelerated particles reveal two kinds of signatures : 
narrow energy spectrum and broad energy spectrum. For instance, Hull et al. [44] show 

an observation where a beam of energetic electrons is seen with a low energy dispersion, 
and a low pitch angle dispersion. This distribution is associated to a localized source of 
acceleration like a strong double layer. Later, a sharp transition occurs. The low energy 
dispersion beam disappears, but we can see energetic electrons with a broader dispersion 
in energy and pitch angle. They are associated to the acceleration by an Alfven wave. 
With double layers, the acceleration region is small, and the accelerated particles arrive 
at the detector with a low dispersion in energy. 

In a case study with FAST data, Chaston et al. [45] analysed an event of electrons 
with a broad spectrum of energy. They where observed in conjunction with Alfven 
waves carrying an important Poynting flux. When particles are accelerated by waves 
of long wavelength such as an Alfven wave, the fastest electrons (they can go faster than 
the wave) quit the acceleration region and run ahead of it. The size of the acceleration 



region is supposed to be of the same order as the wavelength. Because of the time delays 
to reach the detector, and because the acceleration region is extended, the spectrum of 
particles accelerated by a long wavelength wave is broad. A broad spectrum of energetic 
electron is now often used as a proxy for acceleration by a wave. This is the case, for 
instance, in the already mentioned statistical study of Newell et al. [30]. 



The case of yy carried by global field line oscillations 

Linear MHD equations of the oscillations of dipolar magnetic field lines (without 
compressional effect) were derived by Radoski [39]: 

Va X Va X (V X V X E) = y?(V X V X E)^ = -^E. 

dt 

Using dipole coordinates, V = r~^sin^0 and jU = r~^cos0, defining gy = h^E^, e,^ = 
h^E^,H\ = l/h^,H2= l/hy, and E^ = and = 0, for an Alfven wave (perpendic- 

— * 

ular to the dipole field Bq), he found 

V'l^^^^'* = 

where co is the wave frequency, and hy and /i^ are scale factors associated to the dipole 
coordinate system. This system shows poloidal (variations of Cy) and toroidal (variations 
of e,^) vibration of the field lines. Here, they are totally decoupled. These two equations 
have different frequencies for the fundamental mode, but the frequencies are close for 
modes n = 3,.... [40]. Many effects such as plasma compression [41] couple these 
equations and a solution can be found for the fundamental mode and its harmonics 
[42, 43]. Some modes are mostly compressional, others are mostly azimuthal (with no 
compression). Other models include ionosphere's electrodynamic properties [46, 47] 
and finite electron inertia [48]. They reproduce quite well the observations made from 
the ground and from space [49]. Theoretical works show that the irregular Pi2 (<40s) 
could be generated by a surface AW at the inner frontier of the plasma sheet [50] or by 
an AW generated by a sudden layer of anomalous resistivity in the magnetotail [51]. 

Most of the models were developed with MHD or bi-fluid equations. As said before, 
this is not favourable to acceleration studies. Nevertheless, the ability of the FLR to 
accelerate the electrons was investigated by Tikhonchuk and Rankin [52]. As the other 
authors, they based their computation on a MHD model to describe non-compressional 
oscillations of closed field lines. They started from an heuristic equilibrium B field model 
(dipole or tail-like) and ran a MHD numerical simulation, with an initial perturbation of 
the whole field lines. The computation cannot provide a good estimate of but gave a 
parallel current density j , shown in Fig. 3. Then, the authors computed the motion of 
electrons with an initial distribution f{to,v): they bounce between mirror points. Then 




FIGURE 3. Left: Parallel current and mass density perturbation associated to MHD computation of a 
FLR [52]. Right: Parallel electric field in a FLR, according to a purely MHD computation (dotted line) and 
when kinetic aspects of the electron motion are taken into account (solid Une). Reprinted with permission 
from [52]. Copyright [2000], American Institute of Physics. 



they upgraded / at time t, deduced a charge density and and iterated to see 

how, in turn, E\\ accelerate the electrons. The physics involved in this computation is 
that described by Knight [4], Ergun et al. [19], Andersson et al. [25] in the context of a 
direct slowly variable current that is forced by a MHD wave. Tikhonchuk and Rankin 
showed that the reaction of the accelerated electrons forced by the slowly oscillating 
current of the FLR could generate higher electric field than those predicted by a full 
MHD model. Figure 3 (right hand side) shows the parallel electric field provided by their 
MHD model (dotted line), and the larger electric field derived from the computation of 
the electrons motions (solid line). Although it is not totally self consistent, this work 
provided interesting insight on the ability of FLR to accelerate electrons. It is suggested 
that auroras associated to FLR could appear during the growth phase of substorms, when 
the magnetic field lines in the magnetotail are stretched by the intensification of the 
magnetotail transverse current [49]. 

Le Contel et al. [53, 54] developed a kinetic theory of the very low frequency per- 
turbations of a transverse current imposed on a dipolar magnetic field. The perturba- 
tion is supposed to mimic those of the tail current. This current is carried by a plasma 
with ions hotter than electrons. The authors computed the reaction of the plasma, 
with a time varying electric field that is not purely electrostatic 5E = — V54> — dtdX 
but where quasi-neutrality (as in MHD) is relevant. They showed that at zeroth order 
54> + A = ^>o(v^,)^) + 0{Te/Ti), where ^»o(v^,)^) is a function that is constant over a 
magnetic field line, and A is a function related to the vector potential A. Then, at first 
order 5£'|| = — A0) — This electric field generates a FAC system that is mostly 
carried by electrons, dependent on the radial distance, with a sense reversal consistent 
with the observations. 

Full field line resonance implies closed field line. Only early stages of substorms 
(quiet arcs) are concerned. Later, during substorm onsets (first bright and non stationary 
arcs mostly associated to open field lines), smaller scale AW are involved. 



Inertial and kinetic effect that allow for a finite E\ 



An AW can carry a parallel electric field if its perpendicular wavelength is small 
enough, below the limit of validity of MHD. In the low altitude of the auroral zone, 
/3 <^ me/nii. (/3 is the ratio of the plasma pressure to the magnetic pressure.) The motion 
of the particles in the strong magnetic field is given by the guiding-centre approximation, 

-r- = +-E|i and u = Ex-y + — ^— ^. (3) 

Because of their low mass, the electrons have a high mobility along the magnetic field 
lines, while the ions are sensitive to the polarization drift (last term of the second 
equation). When the ions see the perpendicular electric field of the AW, they move 
accordingly across the magnetic field lines. The much lighter electrons cannot do so, 
and a charge density establishes, that causes a parallel electric field which favours 
the resetting of charge neutrality. Therefore, the electrons are accelerated along the 
magnetic field direction. Then [55], the dispersion relation is not exactly those provided 
by the MHD equations, and the parallel electric field becomes important as soon as 
k\c^/(ol >~ 1 

When, 1 /3 ^ nie/mi, the inertial effect are dominated by kinetic effects linked to 
finite temperature [56], 

-^=*fvi(l+^ip?(M)).^ = ^|*ip?. (5) 
A more thorough expansion [57] provides the following dispersion relation 

- >,2x.2 ( kjpf klpKTe/Ti) 



i-r(^ip2) riklp^)[i + i;ziQ] 

where ^ = CO/k^^Ue, Ug — {ITg/mey^^, and Z is the Fried and Conte function. Both 
inertial and kinetic effects can be treated simultaneously, in the fluid approximation with 
a gyroviscous stress tensor [58], leading to 



l+^]_c2/(02 

In all cases, a parallel electric field is set as soon as /3 < mg/m, and k\c^/o)^ ~ 1 

or mg/m, < j8 << 1 and k\pf ~ 1. Figure 1 sketch the regions of the auroral zones 
where these conditions can be met. The inertial region corresponds to inertial AW and 
j8 < me/mu the kinetic region correspond to the regime where me /mi < j8 << 1. The 
plasma j8 = v^j^jV^ is estimated after the profile of the Alfven velocity and the electron 
(or ion) thermal velocity along an auroral magnetic field lines, such as Fig. 4. 




FIGURE 4. Va and Vfhe as functions of the altitude along a magnetic field line Reprinted with permission 
from [43]. Copyright [1995], American Geophysical Union. 



The case of y y carried by mesoscale Alfven waves 

Lysalc and many otlier authors have developed a theory of Alfvenic resonance of a 
field line in the inertial regime k^c/(Op ^ 1, that corresponds to a transverse wavelength 
of typically 10 km in the auroral zone. They assume a non-uniform vertical density 
profile (as in Fig. 4) and they consider that only a portion of the field line is in resonance. 
In that case, the wave frequency is larger than the bounce frequency of the (eventually) 
trapped particles. This simplifies the computations, and it is possible (contrarily to the 
FLR's) to neglect the behaviour of the trapped particles. Thompson and Lysak [60] have 
solved the bi-fluid equations of an inertial AW in an inhomogeneous plasma, in order 
to evaluate the electron acceleration. For a given density distribution (the gradient is 
vertical) they show that a resonant cavity is formed, and the associated frequencies are 
harmonics of a basic frequency close to 2 Hz (see Fig. 5). It is of the order of Va/Lq 
where Lq is the scale height of oxygen. In numerical simulations, they computed the 
electric field, and injected test-electrons, from the upper side of the simulation box. 
They also added a potential drop of 6 kV. This is not self-consistent, but it can help to 
see what kind of electron distribution is produced. They found electron conies similar 
to those observed in the acceleration regions, when the imposed electric field and the 
Alfven wave electric field are of similar amplitude. 

Lysak and Song [61] developed a theory of the resonance of a portion of field line 
when an AW of given frequency (from 0,1 to 1 Hz) and A;^ (in the inertial regime) is 
imposed. This resonance is based notably on the fact that the Alfven velocity has a 

— * 

strong longitudinal (parallel of Bq) gradient at the ionosphere, reaches a maximum, and 
then decrease at altitudes exceeding 1 Ee- This is equivalent to a resonant cavity. It is 
called the ionospheric resonant cavity. They neglected the trapped particles, and this 
theory can be applied to open field lines. They consider that the non-local theory for 
(0> (Oh (where is the trapped electron frequency of bounce between mirror points) 
is relevant for the low altitude part of the magnetic field lines. They show that beyond 
3 Re the inertial effect is negligible, but between 2 and 3 Re the absorption of the wave 



energy by the plasma is important, and this is the place where particle acceleration must 
be expected. At lower altitudes, the wave is reflected. (Simpler models consider the 
wave reflection only at the ionospheric level). A part of the wave energy is dissipated 
in the ionosphere through the Joule heating resulting from the Pedersen conductivity 
(the ionospheric conductivity for horizontal currents that are parallel to the horizontal 
electric field). If the electrons in the magnetospheric side have a temperature exceeding 
a few hundreds eV, a parallel electric field of a few mV/m (of the two signs) can 
be produced. The acceleration produced by these resonant waves happen in both the 
upward and downward directions. Even if auroras are concerned mainly with downward 
acceleration, this does not contradicts the observations: they show fluxes of energetic 
electrons in the two directions. 

The previous works were based on linearised computation, with or without test par- 
ticles. Lysak and Song [62] made a more refined theory, where the electron distribution 
is iterated, and the problem linearised in the vicinity of the parameters of the previous 
iteration. Passing and reflected electrons are taken into account, but trapped electrons are 
still neglected. The results are quite analogous to those of the linear theory (Fig. 5 [62]). 
They show that 10% of the wave Poynting flux can be transferred to the accelerated 
electrons, with an instantaneous rate sometimes reaching 40%. Therefore, the accelera- 
tion might occur as a series of bursts, as was observed with the FAST spacecraft [63]. 
But at this stage, it is necessary to consider the reaction of the ionosphere to the sudden 
input of energetic electrons. (See next section). In this system, the potential drop over a 
wavelength can reach several hundreds eV, therefore, electrons can become trapped in 
the wave. When the wave velocity increases (V^ is not uniform over a field line), trapped 
particles are more efficiently accelerated [64]. 

Clark and Seyler made numerical simulations of inertial AW in the regime 7] >> Tg, 
with a PIC code [65]. Their simulation includes only one wavelength, and consider only 
an homogeneous plasma. There is no convergence of the magnetic field lines, and no 
resonant cavity. But this simulation is self -consistent and fully non-linear. It shows that 
electrons can indeed be trapped in the wave. At the trapping location, the wave profile is 
sharper and trapped particles, initially slower than the Alfven velocity can reach twice 
the wave velocity. 

Watt et al. [66, 67] showed, also with non-linear numerical simulations, that the 

linear conditions of optimal acceleration becomes k±pi > 1 when trapped electrons are 
considered. Therefore, the acceleration is produced by waves of smaller perpendicular 
wavelength than 10 km, when the observations show typical sizes of 10 km. When the 
trapped electrons move into a region where the Alfven velocity is reduced, the electrons 
become de-trapped, and form beams of energetic particles [68]. 

Swift [69] made numerical simulations with a bi-fluid description of the plasma and 
the addition of test-electrons, but in 2D. The convergence of the magnetic field lines is 
taken into account as well as the vertical density gradients and the resonant cavity. The 
wave profile is not uniform over the horizontal axis. He showed that the acceleration 
occurs in the two directions, but is more efficient for upward electrons. He concluded 
that the high energy earthward electrons have been accelerated in the other hemisphere, 
but this supposes that the magnetic field lines are closed. This is also at odds with 
the observations in the upward current region that show that the energetic downward 
electrons are dominant. 




FIGURE 5. Left: Resonance frequencies in an Alfvenic resonator. Reprinted with permission from [60] . 
Copyright [1996], American Geophysical Union. Right: Parallel electric field in an Alfvenic resonator 
caused by an inertial effect with small k±. Reprinted with permission from [62]. Copyright [2005], 
American Geophysical Union. 



Influence of the accelerated electrons on the 
ionosphere-magnetosphere interface 

In the ionosphere, the plasma is coUisional. The electric current Ix flowing horizontally 

is related to the electric field by an Ohm's law, = ^pE^ + "^hEj = N{PEx + HEy) ~ 
NPEx, where Zp and Zfj are the Pedersen and the Hall conductivities, and P and H 
are the corresponding mobilities. The Perdersen conductivity (relating the current to 
the electric field that has the same direction) is the dominant term. The current flowing 
horizontally usually connects regions of downward and upward field aligned currents j. 
The current is divided in two parts, j = j] + j2 where ji is the dominant current carried 
by the downward electrons, and 72 is carried by ions and upward electrons. The flows 
of energetic particles that carry the current contribute to the ionisation of the neutrals, 
modifying the electron density 

AA^ = -j8i7i - i8272 -a{N^- N^) . (6) 

These currents are influenced by the acceleration of downward electrons because of 
their ability, through 72 in Eq. (6) to ionise the neutrals and, consequently, to increase 
the ionospheric conductivity Lp. Furthermore, the current continuity equation is simply 
Ix = J jdx. On the magnetospheric side, the same field aligned current flows, unless a 
part of it goes into a polarization current 



perpendicular to the magnetic field, and induced by the motion of the plasma or by an 
AW. 
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FIGURE 6. Left: Sketch of the relevant parameters for the ionospheric-magnetosphere coupUng. Right: 
Variations of the fields for the ionospheric feedback instabiUty. Reprinted with permission from [70J. 
Copyright [1970], American Geophysical Union. 



The ionospheric feedback instability 

The above equations can be coupled, provided we know the origin of the electric field 
that defines the transverse magnetospheric current. Atkinson has shown that this system 
of current, when it is related to the steady convection of the magnetospheric current can 
generate an instability that increases the field aligned current, and therefore, the ability of 
the plasma to accelerate electrons along the magnetic field [70]. The mechanism of the 
instability is sketched in Fig. 6 (right hand side) and explained in the chain of causalities 
listed below and called the ionospheric feedback instability 

• (a-b) A localised increase of ionospheric conductivity 

• (c) produces an electric field Ex{x) in the ionosphere 

• (d) that is transported into the magnetosphere where it becomes E'^{x). If the 
field line is iso-potential, and if the magnetic field lines convergence is neglected, 

E'^{x)=E,{x). 

• (e) A moving flux tube of velocity Vx (magnetospheric convection) sees E'^ as a 
varying electric field, dfE^ ~ VxdxE'^. 

• (f) It causes horizontal polarization currents defined by Eq. 7 

• (g) closed by field aligned currents j that are added to the initial FAC. 

• (a-b) They intensify the ionospheric ionisation and its conductivity. 

• Goto (a), there is an instability. 



The work of Atkinson has two drawbacks. Atkinson showed that this instability 
generates a weak FAC j when it is only supported by steady convection at stage (e); 



and contrarily to what happens with auroral acceleration, the model does not take into 
account a possible parallel electric field at stage (d) and this is in contradiction with 
our interest for particle acceleration. But this interesting cycle can be improved if we 
consider that the conversion at stage (e) of a spatial derivative of Ex into a temporal 
one is caused by an AW instead of convection, and if the potential drop at stage (d) is 
included in the loop. Trakhtengertz and Feldstein [47] did it, with the consideration of 
the closed field line resonance described by Eq. (2) and with consideration of the work 
of Maltsev on the boundary condition for AW in the ionosphere [46]. They showed that 
the instability cannot occur for a transverse scale smaller than 5 km, and that the motion 
that couples the magnetospheric plasma to the instability (stage (e)) must exceed 100 
m/s. The optimal transverse wavelength is 8km and the corresponding period is about 
10 s. 

This model does not take into account the possible reflection of the wave on the 
ionosphere and thus the ionospheric resonant cavity. The latter was considered by Lysak 
[71]. As it concerns only the low altitude part of the field line, Lysak abandoned the use 
of dipolar coordinates and subsequently, the effects linked to the field line curvature. 
But he included the absorption and the reflection of waves by the ionosphere. For 
the magnetospheric part, he supposed a vertical profile of the Alfven velocity of the 
form = Vai/^s^ + exp—z/h), where z is the altitude, and h is the scale height. 
The electromagnetic field is deduced from the potentials et and the equation of 
propagation of the wave of frequency ft) is 

^ + ft)Vi(z)«J> = 0. (8) 

With the particular profile of {z) given above, and after the introduction of the variable 
X = xoexp{—z/2h), a Bessel equation of negative order is found. The solution takes the 
form ^» = A+7,;i.ye(x) +A^7_jx()e and = {ic /Vai){x / XQ){d^ / dx) . A low order de- 
velopment of the Bessel functions shows that the first term in ^> corresponds to a down- 
going wave, and the other to an upgoing wave. Then 4> = 4>oexp {±.ixQe/2hz — icot). We 
can imagine an ingoing wave and a reflected outgoing wave along an open field line, 
or a closed field line excited from the equator. The ionospheric boundary condition is 
derived from previous works, as 

.■^U+5^*l, = 0. (9) 
In a homogeneous medium, the Alfven wave reflection rate on the ionosphere would be 

R= — = -^—^. (10) 
A+ Ea + Ep 



where Ea = 1 / /ioVk is the wave conductance. With the non homogeneous density profile 
considered by Lysak, this ratio becomes 



where a = Zp/Z^/. At this stage of the computation, Lysak generalised these equations 
taking into account the dependency of the ionospheric conductivity on the FAC j. The 
boundary condition becomes 

/(I ^)d<^dx\,, + '^^\,, = 0, (12) 

xq + iv Lai 

where V depends on the loss of ionospheric electrons by recombination, the ionisation 
is characterised by C7 = {2hQ/VAi){P^±+Hk± x z) -Eq, and is the transverse wave- 
vector of the perturbation and Q a coefficient of proportionality, that depends on Eq. The 
equations are solved numerically. The frequencies and the growth rates are shown on 
Fig. 7, for a realistic range of values of the parameters o and a, for open field lines. 
They allow for instability. 

The source of free energy of the instability comes from the reduction of the dissipation 
of energy in the ionosphere when its conductivity is increased. Lysak [72] has shown that 
a surplus of energy can be injected from the ionosphere into the magnetosphere, and the 
amplitude of the reflected AW can even become larger than those of the incoming AW. 
Then, the amplified wave carries a FAC j that, if the phase relations are favourable, can 
again increase the ionospheric conductivity and the feedback instability. 

Lu et al. [73] have coupled a Hall-MHD model of Field Line Resonance (FLR) [48] 
with a code that couple the auroral precipitations with the ionosphere, the physics having 
been developed in previous works [74, 75, 76]. The Hall-MHD code take inertial effects 
into account, and allows for compressional waves. The magnetic field can be dipolar or 
result from an heuristic magneto spheric model [77]. The authors show that the FLR digs 
auroral plasma cavities, via the ponderomotive force. They are centered at 2500 km of 
altitude, they are 2000 km long, and the plasma density depletion can reach 30%. An 
other simulation done with these codes also reproduced fairly well auroral observations 
of a closed FLR seen by Fast and from the ground. 

Streltsov et Lotko [78] made a series of 2D simulations based on bi-fluid equations, 
with inertial effects, AW reflection on the ionosphere (with the equivalent of the a ratio 
shown above), and the dependency of the ionospheric conductivity on the FAC. An 
arbitrary density profile can be set. They started a series of simulations with an initial 
FAC perturbation on a scale of 120 km and an amplitude of 5 iiA.nT^. In the regime 
>> Ea/, on a time scale of a few minutes, the current perturbation develops fine 
scale structures of 10 km, with amplitudes locally exceeding the initial perturbation by 
an order of magnitude. But these structures do not efficiently propagate at high altitudes. 
In the regime Hp ~ 'Lai, in accordance with the works of Traktenghertz, Lysak etc., 
the small scale structures spread all along the magnetic field line. Practically, it means 
that small scale FAC perturbations can be observed at altitudes of a few thousands of 
km. Upward and downward current regions do not have the same behaviour. Other 
simulations made in the same framework [79] show that upward currents tend to develop 
over a narrower space, but they do not develop the sub-structures shown by Streltsov et 
Lotko [78]. These are the regions that correspond to the bright discrete auroras. These 
regions carry a higher current density than the regions of downward currents. But the 
downward current regions, that are broader, are also those which develop small scale 
structures of typically 10 km. This would be compatible with the fact that the black 



FIGURE 7. Frequency and growth rate of ionospheric feedback instabiHty with AW reflection ( (7 is a 
ionospheric recombination rate, and ^a/^p the conductivity ratio). Reprinted with permission from [81]. 
Copyright [1991], American Geophysical Union. 



auroras, associated to these filamented downward current are narrower tlian tlie bright 
auroral arcs. The role of the high ionospheric conductivity in the development of narrow 
current regions in the region of Alfvenic acceleration has been confirmed in independent 
simulations, when the acceleration was studied in the presence of density gradients that 
are transverse to the ambient magnetic field [80]. We will see in the next section the 
interest of studying the interaction of AW and transverse density gradients. 



Alfven wave + Plasma cavity Electron acceleration + Turbulence 

All the theories of auroral acceleration by AW are based on the existence of small 
transverse scales. What is their origin ? We can suppose that the AW have been generated 
in the magnetosphere, far from the auroral zone, with an already large value of k_i_, and 
transported along the magnetic field lines (the group velocity remaining mostly parallel 
to the magnetic field). It is also possible that the small transverse scale is generated 
locally, precisely at the place where the acceleration occurs. 

The FREJA spacecraft crossed the auroral zone at an altitude of about 1000 km. 
Localized structures called Solitary Kinetic Alfven Wave (SKAW) were observed. They 
consist of isolated electromagnetic perturbation (observed typically during 10 ms) that 



have the characteristic polarisation of AW with small transverse scale (the plasma 
at these altitudes is in the inertial regime j8 < mg/m,) [32]. Some of these structures 
were observed propagating along a plasma cavity with transverse gradients, and in 
the presence of energetic electrons that may have been accelerated by the AW. Their 
Poynting vector is larger than the flux of the energy of the accelerated electrons by 
one or two orders of magnitude [82, 83]; therefore, if one find an acceleration process 
associated to these wave with an energetic efficiency of a few percent, it can be an 
explanation to the observation of energetic electrons in the vicinity of the SKAW's. 

Let us consider AW initially set with = 0, propagating along the magnetic field 
direction. Initially a wave front is a horizontal surface. But in an auroral cavity, the 
magnetic field is hardly modified when the density can reach a very deep depletion, 
as was shown experimentaly [84] and theoretically [85]. The transverse gradients can 
be very sharp, with a typical size comparable to a Larmor radius that, in this area, is 
also comparable to the electron inertial length. Therefore, the Alfven velocity Va — 
B/ {iiqPm) where Pm is the mass density, increases in the low plasma density region. 
Then, the wave propagates faster in the cavity than outside, and the wave front is 
distorted. It becomes oblique to the magnetic field along the transverse density gradients, 
with the development of small scales comparable to those of the transverse gradients, 
that are of the order of the inertial length [86]. As was shown with Eq. (4), when k± 
reaches this scale, the AW develops a parallel electric field, and acceleration becomes 
possible. Of course, the computation of (4) is based on a homogeneous plasma. Genot et 
al. [87] derived the parallel electric field of an AW propagating upon a density gradient. 
Their computation is based on the linearised bi-fluid equations of a cold plasma. They 
include shear AW and inertial effects. They found, when the gradient is purely transverse 
that, 

Ell _ (dxn/nfc^/(oj 
E± ~ (dxn/n) 1 + {dxn/n)^c^/(oj ' 

This relation is analogous to Eq. (4), provided k± is replaced by dxn/n. 

But the work of Genot et al [87] is based on linearised equations, and does not 
describes the kinetic behaviour of the electrons. It is therefore inappropriate for studying 
the particle acceleration associated to this electric field, and cannot show the retro-action 
of the accelerated electrons on the plasma that proved so important in the works cited up 
to now. Therefore, a series of numerical simulations was conducted with a PIC code [88] , 
to study the propagation of an AW upon a purely perpendicular sharp density gradient. 
The simulations were initialised with a deep plasma cavity (generally dn/n = 0.75) 
with gradients developing over a few electron inertial lengths, and the superposition of a 
shear or magnetosonic AW propagating along the ambient magnetic field. As predicted 
by Genot et al. [87], parallel electric fields with the same wavelength as the AW develop 
along the density gradients. These electric fields, in the regions where they are strong, 
heat the electrons and some of them reach velocities greater than Va- Then, the most 
energetic electrons escape from the high parallel electric field regions, forming beams 
[89, 90]. The electron beams quickly generate charge separation (the ions do not follow 
the electrons) that result in electrostatic instabilities. A complete re-organization of the 
electric field follows the growth of the electrostatic instabilities, with localized electric 



structures of high intensity (locally much higher than the AW parallel electric field of 
Eq. (13)). These localized spikes of electric field are associated to vortices in the electron 
phase space projected along the magnetic field direction (x||,V||). These structures were 
identified as weak double layers and electric solitary waves [91] . The acceleration was 
associated to the deformation of the Aflven wave front and to the creation of small scales. 
The turbulence generated by the accelerated electrons is analogous to that observed 
in the auroral zone in the vicinity of the auroral plasma cavities. It was shown that 
when this turbulence is fully developed, the acceleration by the AW looses its efficiency. 
Acceleration by AW starts again when the weak double layers have dissipated [92]. 
These studies were conducted first with a sinusoidal AW. But the small scale AW 
observed in the auroral zone are not monochromatic, they are SKAW's. Therefore, a 
series of numerical simulations was conducted with an Alfven wave packet. It was 
shown that the results found with a monochromatic wave survive with isolated wave 
packets. It was also possible to show that the AW acceleration increases by an order of 
magnitude the flux of kinetic energy of the accelerated electrons, that is the parameter 
that characterize the ability of the energetic electrons to trigger the auroral visual display 
(excitation of oxygen and nitrogen atoms in the ionosphere)[93]. 

The interaction of AW with a perpendicular density gradient were studied with SKAW 
because, apart from being observed, they can be simulated with a PIC code thanks to 
their relatively short wavelength. But the interaction of AW with density cavities can 
also happen on a larger scale. Lysak and Song [80] made 3D simulations of mesoscale 
AW trapped in the ionospheric Alfvenic resonator when it includes a plasma density 
depletion (and subsequently transverse density gradients). This study was conducted 
with a bi-fluid code and it cannot describe the electron acceleration and the subsequent 
plasma turbulence generated by the energetic electrons. It nevertheless showed that the 
AW develop small transverse scales, and that parallel electric fields caused by inertial 
effects grow along the density gradients. 

DIVERSITY OF AURORAL ACCELERATION PROCESSES 

This review has shown that the auroral acceleration can be caused by: (i) Quasi-static 
electrostatic structures called double layers. They are a combination of the field line 
convergence effect, and of the retro- action of the accelerated particles on the plasma that 
tends to localise the electric potential jump over a reduced distance along the magnetic 
field direction; (ii) Magnetic field line resonance. (Large scale Alfven waves); (iii) 
Alfven waves trapped in a resonator, i.e. a region around the locus of the maximum value 
of Va- Meso-scale; (iv) Isolated and small Alfven wave packets, that become sources of 
plasma acceleration when they propagate through transverse plasma density gradients. 

The MHD waves do not carry a parallel electric field, therefore non-MHD effect are 
required for electron acceleration. They are (i) Electron inertial effects when k±c/(Ope ~ 
1. They are efficient for acceleration at lower altitudes, where j8 < me/mf, (ii) Ion 
temperature effect when kj^Pi ~ 1 are efficient at higher altitudes where the plasma beta 
me /mi <C j8 < 1. In both cases, this means a need for small transverse scales. They can 
be created by (i) a source of Alfven waves, especially for FLR; (ii) by cCoupling of the 
waves with the ionosphere, in the ionospheric feed-back instability, especially efficient 



for meso-scale AW trapped in the ionospheric AW resonator; (iii) by propagation upon 
a density gradient. Their efficiency was proved for small scales AW with PIC numerical 
simulations. 

New topics are emerging concerning auroral plasma acceleration. For Earth auroral 
processes, many recent works are devoted to the relation between the quasi-static accel- 
eration structures and the acceleration by AW. Up to now, most of them are based on 
observations. Some observations tend to show that the ion beams that have been accel- 
erated by strong double layers can trigger AW instabilities that in turn could generate 
Alfvenic electron acceleration [94, 95, 96]. On the contrary, other observations tend to 
show that Alfvenic acceleration acts as a precursor of the acceleration by strong double 
layers[97, 98, 44]. 

Auroral acceleration processes are also investigated on other planets. Observations 
are generally no done in-situ, because the few orbiters around the other magnetised 
planets (Galileo, Cassini) do not go across the auroral zone. The works on auroral 
acceleration are mainly supported by observation of auroral display (from IR to UV) 
[99, 100] and of the associated radio emissions. For Jupiter, the radio emissions coming 
from auroral regions are in the decametric range and can be observed with ground 
based radio-telescopes. The interpretation of decametric radio emission from the flux 
tube that connects Jupiter to its closest big satellite lo has been interpreted in term of 
sources emitting at the local electron gyrofrequency by bunches of energetic electrons 
[101, 102]. It has been shown that the quasi-periodic character of these emissions and 
other characteristics of the dynamic spectrum of these emissions can be due to the 
presence of quasi-monochromatic AW in the auroral region [103, 104, 105] that would 
accelerate downgoing electrons [106, 107, 108, 109]. The in-situ exploration of the 
auroral zone of Jupiter is now in project [110], and theoretical models already allow 
to speculate on what a probe would observe in the acceleration regions [1 1 1]. 
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